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ABSTRACT 

Aims. This work aims to provide a first insight into the mass-metallicity (MZ) relation of star-forming galaxies at redshift z ~ 1.4. 
To reach this goal, we present a first set of nine VVDS (VIMOS VUT Deep Survey) galaxies observed with the NIR integral-field 
spectrograph SINFONI on the VLT. 

Methods. Oxygen abundances are derived from empirical indicators based on the ratio between strong nebular emission-lines (Ho-, 
[Nii]6584A and [S ii]67 17,673 lAA). Stellar masses are deduced from SED fitting with Chariot & Bruzual (2003) population synthesis 
models, and star formation rates are derived from [Oii]3727A and Ha emission-line luminosities. 

Results. We find a typical shift of 0.2 - 0.4 dex towards lower metallicities for the z ~ 1.4 galaxies, compared to the MZ-relation in 
the local universe as derived from SDSS data. However, this small sample does not show any clear correlation between stellar mass 
and metallicity, unlike other larger samples at different redshift (z ~ 0, z ~ 0.7, and z ~ 2). Indeed, our galaxies lie just under the 
relation at z ~ 2, and show a small trend for more massive galaxies to be more metallic (~ O.I logarithmic slope). There are two 
possible explanations to account for these observations. First, our galaxies present higher specific star formation rates when compared 
to the global VVDS sample which could explain the particularly low metallicity of these galaxies as already shown in the SDSS 
sample. Second, inflow of metal-poor gas due to tidal interactions could also explain the low metallicity of these galaxies as three of 
these nine galaxies show clear signatures of merging in their velocity fields. Finally, we find that the metallicity of four galaxies is 
lower by ~ 0.2 to 0.4 dex if we take into account the N/O abundance ratio in their metallicity estimate. 

Key words, galaxies: abundances - galaxies: evolution - galaxies: high-redshift 



^ 1. Introduction iGallazzi et all 120051: [Paiiter et all 12008*) and the 2 degree Field 



Galaxy Redshift Survey (2dFGRS; Lamareille et al. 2004). The 



The abundance of heavy elements in galaxies reflects the past relation for local galaxies is steep for masses < 10»'^ 

history of star formation and the eff^ects of gas exchanges (in- g^^^^^ ^^^^^^ ^^^^^^ suggesting that the MZ rela- 

flows and/or outflows) with the intergalactic medium. A charac- jj^^ ^^^j^jy ^^j^g^ decrease of metal loss when stel- 

terization of the evolution of chemical abundances for galaxies ^^^^ increases. These trends observed by iTremonti et all 

of different masses is therefore essential to a complete model ^2004) and Lamareille et_al. (2004) have been shown to extend 

of galaxy formation that includes the physics of baryons (see m uch lower galaxy m asses (M < 10'' Mo; eg. Lee et al. 

egjDe Lucia et al. 2004; de Rossi et al. 2007; Finlator &_Dave 2006: ISaviane et dlM . Recent studies focused on the depen- 



22221' ■ dence of the MZ relation of SDSS galaxies with environmen t 

First discovered for irregular galaxies (| Lequeux et al. || 1979D , „„ g^ales ("Michel-Dansac et al. 2008: ElHson et al. 2008bh 



on small scales (Michel-Dansac et al. 2008; Ellison et al. 2008b) 
the mass-metallicity (hereafter MZ) relation has been investi- j^gg^ g^^jg^ ( Mouhcine et al. 2007; Cooper et al. 20081) 



gated intensively (Brodie & Huchra 1981; Skillman^etalJll989t and th eir star formation rate and galaxy sizes (lElhson et alJ 

iZaritskv et al.i.l994; Garnett et al.lll997l among others) and is |2008al) 
now well estabUshed in the local universe iz < 0.2) thanks 

to the analysis of data from large spectrosco pic survey s such The evolution of the MZ relation on cosmologicaltimescales 

as the Sloan Digital Sky Survey (SDSS; Tre monti et al]|2504t is now predicted by semi-analytical models of galaxy formation. 

Reliable observational estimates of the MZ relation of galax- 

* Based on observations collected at the European Southern ies at diff^erent epochs (and hence diff'erent redshifts) may thus 

Observatory (ESO) Very Large Telescope, Paranal, Chile, as part of the provide important constraints on galaxy evolution scenarios. 

Programs 75. A-03 18, 78. A-0 177, and 070. A-9007 At increasing redshifts, as the strong rest-frame optical emis- 
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sion lines shift into the near-infrared (NIR) window, metallic- 
ities are typically based on smaller subsets of strong emission 
lines, through the use of empirically calibrated abundance in- 
dicators. Much progress has been made recently in assembling 
large samples of star-forming galaxies with abundance measure- 
ments at both intermediate redshifts (z < 1; Savaglio et al. 2005; 
Maier et al.' '2005'; 'Liang et al.' 2006; Lamareille et al. i2006i 
2009; Perez-Montero et al. 2009; Cowie & Barger 200 83 and at 
Z>2 (Erb et al. 2006; MaioHno et al. 2008; Ha vashi et al.l2009l: 
iMannucci et al..,2009.) . However, chemical abundance measure- 
ments ar e only available for a limited number of galaxies a t 
z~ 1 -2 (IShaplev et al .120051; iMaier et al.l2006tlLiu et al.l2008h . 
The rather unexplored 1 < z < 2 redshift regime is one of partic- 
ular importance in the history of the universe corresponding to i) 
the peak of star formation rate, and hence metal product ion, for 
the universe as a whole (see eg. Lilly et al. 1996; Mad au et al.l 
119961; iPerez-Gonzalez et al.l 12^51; ITresse et alj 120071) . ii) the 
buildup o f a significant fraction of the stellar mass in the uni- 
verse (eg. iDrorv et al.1 120051: lArnouts et all 120071: iPozzetti et all 
120071) . and iii) the emergence of the Hubble sequence of disks 
and spheroids. 

We have undertaken an large observing program at ESO- 
VLT aimed at probing the mass assembly and metallicity evo- 
lution of a representative sample of galaxies at 1 < z < 2. 
The final goal of the MASSIvQ (Mass Assembly Survey with 
SINFONI in VVDS) project is to obtain a detailed description 
of the mix of dynamical types (rotating disks, spheroids, and 
mergers) at this epoch and to follow the evolution of fundamen- 
tal scaling relations, such as the Tully-Fisher and MZ ones, and 
therefore constrain galaxy evolution scenarios. In this paper we 
present the first results obtained by such a program focusing on 
the M Z relation at z ~ 1 .4. A companion paper (lEpinat et al.l 
|2009|) is devoted to the kinematical analysis of galaxies at these 
redshifts using the same data. This work is based on a sam- 
ple of n ine galaxie s selected in the VIMOS VLT Deep Survey 
(VVDS; lLe"Fevre et al. 2005) and observed in the NIR with the 
3D-sp ectrograph SINFONI on the ESOA'LT (lEisenhauer et alj 
l2003h during pilot observing runs. Throughout this paper we as- 
sume a standard A-CDM cosmology with h = 0.7, Qa = 0.7, 
and Q,„ - 0.3. 



2. Data description 

2.1. Observations 

The selection of t he sample and the S INFONI observations are 
fully described in lEpinat et alj (l2009t) . For the convenience of 
the reader, we summarize here the main steps. 

We have used the VIMOS VLT Deep Survey to select galax- 
ies with know spectroscopic redshifts in the range 1 < z < 2. The 
VVDS is a complete magnitude selected sample avoiding the bi- 
ases linked to more crude a priori color selection techniques. It 
off'ers the advantage to combine a large sample with a robust se- 
lection function and secure spectroscopic redshifts necessary to 
engage into long single objects integrations with SINFONI be- 
ing sure to observe the Ha line outside of the OH sky emission 
lines. 

The nine galaxies studied in this paper have been selected 
among those showing the strongest [On] 3727 A emission line 
(EW > 50 A and flux > 5 X 10~'^ ergs^'cm"^, as measured 
on VIMOS spectra) for Ha to be easily detected in the near- 
IR. These criteria for selecting late-type star-forming galaxies 



have been shown to be very efficient. From the first observ- 
ing runs, our success rate of selection has been around 85%: 9 
galaxies over 12 observed show strong rest-frame optical emis- 
sion lines in SINFONI datacubes. The other three galaxies have 
not been detected with SINFONI for the following reasons. The 
redshift of one galaxy is such that Ha is falling just in the gap 
between the J and H bands. For another galaxy, the strength of 
[O ii]3727A was overestimated. New measurements of this emis- 
sion line show that this galaxy is now outside of the selection 
box. The third non-detected galaxy had a wrong redshift. 

Among the nine VVDS star-forming galaxies studied in this 
paper, five have been selected in the VVDS-22h wide field 
(17.5 < Iab < 22.5) and four in the VVDS-02h deep field 
(17.5 < Iab < 24.0). Their basic properties (RA, DEC, red- 
shift, /-band magnitude, etc) are listed in Table [T] These targets 
span a redshift range between 1.27 < z < 1.53. The NIR spec- 
troscopic observations were acquired with the 3D-spectrograph 
SINFONI at ESO-VLT during two 4-nights run on September 5- 
8, 2005 (ESO run 75.A-0318) and November 12-15, 2006 (ESO 
run 78.A-0177). SINFONI was used in its seeing-limited mode, 
with the 0.25" x 0.125" pixel scales leading to a field-of-view 
of 8" X 8", and the H-band grism providing a spectral resolu- 
tion R ~ 4000. Note also that two galaxies (VVDS220544103 
and VVDS220584167) have also been observed with the J- 
band grism. Conditions were not photometric and the mean see- 
ing was around 0.65". The total on-source integration times, 
as well as general information are listed in Table [T] For each 
Observation Block a hot telluric standard star has been observed. 
With this star we were able to flux-calibrate our observations 
with an accuracy of about 20%. Details for data reduction are 
available in ,Epinat et aL (,2009) . 

2.2. Measurement of emission lines 

To perform the spectrophotometric analysis of the reduced dat- 
acubes and measure the detected emission lines, we first pro- 
duced ID spectra from the 3D data. It has been performed in 
two steps: first, we produced a 2D pseudo-slit spectrum from the 
datacube (as if we placed a long-slit over the spatial extent of 
the object), and then we reduced the 2D-spectrum in a ID-one, 
taking care of the spatial extent of the galaxy and of the evolu- 
tion of its position through the dispersion axis. The first step has 
been done with the QFitsView softwar^ which is a tool to vi- 
sualize and analyse datacubes. The second step was performed 
using the IRAF software, specifically applying the apall task 
of the noaotwodspeoapextract package. Figure [1] shows 
examples of ID spectra for the galaxies VVDS020147106 and 
VVDS220015726. 

At the distance of the sample galaxies, the Ha, [Nii]6584A 
and [Sii]6717,6731A fines are redshifted into the SINFONI H 
band. We were thus able to measure these lines in quasi-all our 
sample galaxies. For two galaxies we had observations in the 
J band, which enables to detect and measure [O iii]4959,5007A 
and H/3 lines for VVDS220544103. Unfortunately, the redshift 
of the second galaxy VVDS220584167 makes the reliable mea- 
surement of these lines impossible as they fall at the location of 
the bright OH-lines. We measured the flux in the different lines 
fitting a gaussian to the ID spectrum of each galaxies (within 
the splot task in IRAF). Fluxes and errors are obtained with the 
deblending mode using 100 simulations for the poissonian noise 
model (see the splot task help for more information). Results 
are shown in Table |2l Two reasons can account for the impos- 
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Table 1. Properties of the nine galaxies observed with SINFONI. 
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Galaxy 


RA(2000) 


DEC(2000) 


fexp (band)' 
[hours] 




ZviMOS*^ 


ZsiNFONI^ 


VVDS020116027 


02:25:51.2 


-04:45:07.5 


1.7(H) 


22.87 


1.5259 


1.5302 


VVDS020182331 


02:26:44.3 


-04:35:51.9 


3(H) 


22.73 


1.2286 


1.2283 


VVDS020147106 


02:26:45.5 


-04:40:46.9 


2(H) 


22.50 


1.5174 


1.5195 


VVDS020261328 


02:27:11.1 


-04:25:30.7 


1(H) 


23.90 


1.5291 


1.5289 


VVDS220596913 


22:14:29.2 


+00:22:15.7 


1.75(H) 


21.84 


1.2667 


1.2662 


VVDS220584167 


22:15:22.8 


+00:18:47.5 


1. 75(H), 1 (J) 


22.04 


1.4637 


1.4659 


VVDS220544103 


22:15:25.5 


+00:06:40.0 


1(H),1(J) 


22.47 


1.3970 


1.3966 


VVDS2200 15726 


22:15:42.5 


+00:29:04.9 


2(H) 


22.47 


1.3091 


1.2930 


VVDS220014252 


22:17:45.5 


+00:28:40.3 


2(H) 


22.10 


1.3097 


1.3101 



' Exposure time corresponding to the effective time spent on the source in each band (H or J). 

^ Redshift based on the VIMOS spectrum. 

' Redshift estimate from the position of Ha in SINFONI spectra. 
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Fig. 1. Examples of ID spectrum for two galaxies 
VVDS220015726 (bottom) and VVDS020147106 (top) 
observed with SINFONI in the H band. The spectral re- 
gion of interest is centered around the Ha, [Nii]6584A and 
[S ii]6717,6731A emission lines. The position of these lines is 
indicated by vertical dotted lines. The other lines with P-cygni 
profiles visible in the spectra are residuals of OH sky lines 
subtraction, they are shown within grey zones 



sibility to measure a line: either an OH sky line falls right upon 
the line of interest, or the intensity of the line is too low to be 
detected. Over the nine galaxies, [Nii]6584A is not detected for 
two (VVDS020261328 and VVDS020182331). 

The measured values of log([Nn]6584/Hc!') range from -1.3 
to -0.64, which indicates that the sources are star-forming galax- 
ies with at most a low AGN contamination. We also checked for 
possible X-ray identification of AGN in our sample using the ex- 
tensive HESEARC archive. Neither XMM nor Chandra source 
has been identified in the vincinity of our targets further exclud- 



ing any AG N contaminat i on in our sample. This issue is further 
discussed in lEpinat et al.l (l2009h . 



3. Stellar mass and metallicity estimates 

3.1. Stellar mass 

The procedure to derive stellar mass of our galaxies is based on 
the comparison between observations, basically the photometric 
SED of galaxies, and population synthesis models. More pre- 
cisely, M* estimates were obtained with the G OSSIP spectral 
energy distribution (SED) modeling software dFranzetti et al.l 
[2008). For the SED fitting the multi-band photometric observa- 
tions available in the VVDS fields were used, including BVRI 
data from the CFHT, UBVRZs data from the CFHT Legacy 
Survey, J- and A'-bands data from SOFI at the NTT and from the 
UKIDSS survey, and the VVDS-Deep spectra. The stellar popu- 
lation models used for the fit were generated with the BC03 pop- 
ulation synthesis code (Bruzual & Chariot 2003), assuming a set 
of "delayed" star formation histories (see iGavazzi et a n l2002L 
for details), and a Salpeter IMF. The values and erros bars of 
the stellar masses in Table [3] are the median of the Probability 
Distribution Function and its confidence regions. More details 
can be found in lEpinat et al.l (i2009h . 



3.2. Metallicity 

To determine the gas-phase oxygen abundance of the nine 
VVDS galaxies, we made use of well-known empirical abun- 
dance indicators. Indeed, at these redshifts, we only have access 
to the brightest recombination and forbidden emission lines of 
the ionized interstellar gas, preventing any accurate metallicity 
estimate based on faint [Om] auroral lines. Different calibra- 
tions or empirical indicators have been proposed over the past 
20 years, each of them taking into account different set of emis- 
sion lines. The most commonly used, but which has the draw- 
back of being degenerated, is certainly R23 = ([On] 3727 + 
[Oiii]4959,5007)/H6, first introduced by Pagel et al. (1979), 
and re-calibrated later in many art icles (Kewley & Dopita 2004: 
iNagao et al.ll2006tlYin et al.1l2006h . 

Considering that the nine galaxies have been observed with 
SINFONI in the H band, we had access to Ho- and [Nii]6584A 
emission-line fluxes for six galaxies (we did not detect [N ii]6584 
for three galaxies). This led us to use the N2 indicator, defined 
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Table 2. Emission-line flux (in s"' cm"^) for the nine galaxies observed with SINFONI. 



Galaxy [Oii]3727' H/3 [Oiii]4959 [Oiii]5007 Ha [Nii]6584 [Sii]6717 [Sii]6731 



VVDS020116027 


6.03 ± 


1.11 






6± 


1 


1 ± 


0.3 






0.7 ± 0.4 


VVDS020182331 


19.85: 


t 1.41 






55 ± 


14 












VVDS020147106 


8.53 ± 


1.52 






22: 


t6 


2.9 


± 1 


1.4 ± 


0.6 




VVDS020261328 


10.20: 


t 2.01 






5± 


1 












VVDS220596913 


18.55: 


t 1.78 






33: 


t6 


4.7 


± 1 








VVDS220584167 


17.90: 


t 1.83 






48 ± 


10 


11 


+ 3 


6± 


3 


7±3 


VVDS220544103 


25.32: 


t 2.97 38 ± 9 


46± 12 


114 + 14 


161 : 


t29 












VVDS2200 15726 


9.47 ± 


1.41 






26: 


t3 


4.6 


± 1 








VVDS220014252 


37.54: 


t 1.73 






212: 


t27 


48 


±6 


15 H 


: 2 





' as measured in the optical VIMOS spectra 



as follows: 



N2 = log 



[Nn]6584 



(1) 



which is subject to many calibrations for metalli c ity estimates 



jKewlev & Dopit al 120021; iDenicolo et all l20a 



. - 2t Nagao 

120061; I Yin et alji2006l; iPerez-Montero & Continil 12009) In our 
cases, the N2 parameter offered two main advantages to 
derive metallicity. Due to the proximity of the Ha and 
[Nn]6584A hnes, the N2 ratio is independant of reddening cor- 
rection and uncertainty in relat ive flux calibration. Recently, 
iPerez-Montero & Continil (|2009|) (hereafter PMC08) provided a 
new calibration of this parameter based on a sample of objects 
with a direct determination of the corresponding ionic abun- 
dances: 

12 H- log 0/// = 9.07 -H 0.79 xN2 (2) 

Nevertheless, this parameter presents a large scatter, 
which can be due to man y processes, as explained in 
iPerez-Montero & Diazl ( l2005h . One possible explanation is 
the dependence of the calibration on the nitrogen-to-oxygen 
abundance ratio N/0. PMC08 investigated how the N2 abun- 
dance indicator can be altered by the value of the N/0 ratio. 
They found the following deviation: 



et al.1 



A(0/H) = 0.5 log A^/O + 0.66 



(3) 



Following their method, we had to determine the r iitrog en-to- 
oxygen abundance ratio. 'P erez-Montero & Continil ( 1200 9) pro- 
pose an empirical calibration for it, derived from the same sam- 
ple of H n regions they used to derive equations dU and (O. For 
our purpose, considering the emission lines we have in our spec- 

h H ti, 1 [Nii]6584 
tra, we have used use the log 

[Sii]6717,6731 
ing to the following cahbration: 



N2S2 ratio lead- 



\og(N/0) = 1.26 X N2S2 - 0.86 



(4) 



In this work, we first calculated the metallicities using equa- 
tion (|2]l: the results are listed in Table[3] Unfortunately, the mea- 
surement of the sulfur lines ([S ii]67 17,6731 A) is very difficult, 
and we sometimes only had a measurement for one of the two 
lines, so we derived the other using the relation Q which is valid 
if we assume a typical density of 100 cm"^ (iRamsbottom et all 
Il993h : 

t^^ = 1.32 (5) 
[Sii]6731 

In lOsterbrocld ([T97l . Fig 5.3 p. 134, we can see that at these 
typical densities, the value of the ratio changes by 0. 1 at most. 
We were thus able to "correct" the metallicities of 4 galaxies, the 



ones for which we had a measurement of [S ii]6717,6731A. The 
derived metallicities, corrected for the N/ O ratio, are also listed 
in Tabled 

For three galaxies (VVDS220544103, VVDS020261328 
and VVDS020182331) we did not measure the [Nii]6584A 
line. However, observations in the J band were available 
for the first galaxy (VVDS220544103) enabling to measure 
[Oiii]4959,5009A and HyS emission Hnes (see Table EJ. With 
Ho-, measured in the H band, and [O ii]3727A (from the VIMOS 
spectrum) we had thus in hand five emission lines for this 
galaxy to derive metallicity using the /?23 indicator The ratio 
of the Balmer lines Ha/H/? gives an estimate of the attenua- 
tion by the dust. Assuming a classical value for the intrinsic 
ratio F(Ha)/F(lip) of 2.85, our measurements of the Ha and 
H/5 lines gave an extinction coefficient of C = 0.31 + 0.33. We 
were then able to correct the emis sion-line fluxe s for reddening 
using the Calzetti extinction law (ICalzettill2001i) and use these 
values to derive metallicity. We obtained log 7^23 = 0.702 which 
led to an estimate of the metalli city of 12 + \og{0/ H) — 8.78, 
using the calibration of Zaritsky et al.l (1199 4). Kewle v~& EUisonl 
(l2008b provide a way to transla t e thi s value to T04 (the cal- 
ibration used in iTremonti et al.l (12004.) ). We finally obtained 
12 + log(O/i/)|T04 = 8.74 for VVDS220544103. 



4. The Mass-Metallicity relation 

4. 1 . The MZ relation 

Figure |2] shows the position of the seven galaxies at z ~ 1 .4 in 
the MZ-plane (blue circles) using the metallicities and the stel- 
lar masses listed in Table [3] In order to compare with the MZ 
relations previously derived at different redshifts, we also draw 
in this figure the MZ relations obtained in the local Universe 
(Tremonti et al. 2004) as derived fro m the die SPSS, at z ~ 0.7 
dLamareiUe et al.ii2009t) and z ~ 1 (iPerez-Montero et al.ll2009l) 
from die VVDS sample, at z ~ 1 from zCOSMO S (Contini 
et al., in prep), and at z ~ 2 dErb et al.l |2006). However, 
su ch a comparison r equires great precautio ns, as pointed out 



Kewlev & EUisonl (|i008). In their study, iKewlev & ElhsonI 



(120081) " investigate the influence of different metallicity calibra 
tions on the MZ relation and showed that the absolute value of 
the metallicity can vary up to 0.8 dex depending on the cali- 
bration used. Taking this remark into consideration, our purpose 
being to compare the MZ relation at z ~ 1 .4 to the ones at oth- 
ers epochs, we had to be consistent with the calibrations used in 

other studies. 

Firstly, at z ~ iTremonti et al.l ( |2004|) used an 
indicator based on the Bayesian calibration CLOl (for 
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Table 3. Metallicity (12 + log{0/H)) and stellar mass estimates for the 7 galaxies observed with SINFONI. 



5 



Galaxy 


ZpMCOS ' 


ZpMCOScor.^ 


ZpMCOS-.™' 


ZpMC08cor.^T04^ 




logWO)^ 


log(SSFRyr-')'5 


VVDS020116027 


8.46 ± 0.09 


8.43 ± 0.27 


8.66 ±0.12 


8.62 ±0.37 


io.o9+;j-^» 


-1.13 


±0.49 


-8.72 ± 0.20 


VVDS020147106 


8.37 ±0.16 


8.15 ±0.16 


8.56 ±0.2 


8.29 ±0.19 


10 23+"-'" 


-0.77 


±0.22 


-8.22 ±0.31 


VVDS220596913 


8.40 ±0.11 




8.59 ±0.14 




10.93ii 






-8.94 ±0.21 


VVDS220584167 


8.56 ±0.13 


8.44 ± 0.24 


8.80 ±0.18 


8.64 ±0.33 


11.08!°:^ 


-0.95 


±0.43 


-8.73 ± 0.24 


VVDS220544103 










1Q 71+0.16 






-7.81 ±0.21 


VVDS2200 15726 


8.44 ±0.11 




8.63 ±0.13 




10.79+»-f9 






-8.88 ±0.13 


VVDS220014252 


8.56 ±0.08 


8.20 ±0.12 


8.80 ±0.10 


8.35 ±0.15 


10.78!o;«^ 


-0.53 


±0.22 


-7.83 ± 0.30 



' metallicity estimate using PMC08 indicator. 

- metallicity estimate using PMC08 indicator and N/0 correction with the sulfur lines. 
^ corresponding metallicities converted into T04 metallicity estimator. 
■* stellar mass in solar units (see § 13. It . 
^ log(A'/0) estimated with the N2S2 ratio. 
^ specific SFR in yr"', derived with L(Hq'). 



ICharlot & Longhettil2001h - we will refer as T04 - which make 
use of a large set of emission lines to derive metallicity. In or- 
der to transl ate our metallicities int o the T04 calibrator, we first 
converted to iDenicolo et alj (|2002|) N2 calibration and then ap- 
plied the transformation given in iKewlev & EllisonI (|2()08') to 
get T04 metallicities. If x stands for the metallicity estimate of 
iPerez-Montero & Continil ( |2009|) and y for that of T04, we then 
have: 



-0.217469x^ + 5.80493x2- 



50.2318X-H 149.836 (6) 



The resulting points in Fi g.[2]are plo t ted as blue circles. 

Secondly, at z ~ 2 lErb et aTl ( 20061) used another cali- 
bration of the N2 indicator: PP04 (" Pettini & Pa"gef 2004). To 
be consistent we c hose to convert their relation into T04. In 
IKewlev & EUlsoiil (l2008l) we found the way to convert from 
PP04 to T04, as follows (if x stands for the estimate of metal- 
Hcity with PP04 and y with T04): 

y = 2.6766690x^-68.471750x^ + 585. 1750x- 1661.9380 (7) 

The result ing curve is shown in F i g. [2] in red. The three other 
relations (iLamareille et all 120091; iPerez-Montero et al.1 l2009l 
Contini et al. in prep.) are already caUbrated with T04 estima- 
tor. 

While 'Tremonti et all (|2004. 'Erb et alj (l2006h . 



ILamareille eta l. (2009) and iPerez-Montero et al. (2009J) found 
a monotonic relation between the metallicity and stellar mass 
- more massive galaxies have a higher nebular metallicity than 
less massive ones - our small sample of galaxies at z ~ 1 .4 does 
not show any obvious correlation between metallicity and stellar 
mass, even though one could imagine Z to increase with stellar 
mass. Actually, if we fit a line through our 7 points, we obtain 
a positive slope (12 + log{0/H) = 0.101 log(M*/M0) + 7.61). 
The mean metallicity for this sample of seven galaxies is 
12 + log(0/H) = 8.62 with a scatter of 0.22 dex. Note however 
that the z ~ 1 .4 galaxies are lying below the MZ relation derived 
in the local universe from SDSS data (blue line in Fig. |2]i. The 
shift in metallicity is roughly -0.3 dex, which clearly confirms 
the evolutionary trend of the galaxies metallicity with cosmic 
time for a given stellar mass. Six out of seven gala xies have 
metall icities lower than the median z ~ 2 relation of l Erb et"al1 
(I2OO6I) which seems to draw an upper limit in metallicity for our 
sample. We can wonder why our galaxies seem poorly metallic, 
and in particul ar less metal lic than Erb et al. (2006) relation. 

Recently iElhson et al.l (l2008ah investigated the systematic 
effects of specific star formation rate (SSFR) on the MZ relation. 




log(MVMe) 

Fig. 2. Mass - metallicity relation for the sample of seven 
VVDS galaxies (blue cuxles) at z ~ 1.4. The arrows show 
how the metallicities change when we include t he co rrection of 
log(N/0), following IPerez-Montero & Continil (12009 ) prescrip- 
tion. The error bars at the bottom right corner show typical un- 
certainties of our data. The blue curve is the MZ relation in 



the local universe as derived from the SDSS by .Tremonti et al 
(2004). The red curve is the MZ re lation at z ~ 2 bylErb etal 
(i2QQd) . The MZ relation at z ~ 7 (ILamareille et al.1l2009t) and 
z ~ 1 ( iPerez-Montero et al.ll2009h from the VVDS sample, and 
at z ~ 1 from zCOSMOS (Contini et al., in prep.) are shown in 
black open triangles, black squares, and magenta circles respec- 
tively. All metallicities have been translated into T04 to facilitate 
comparisons. 



At a given mass, SDSS galaxies with high SSFR have in aver- 
age lower metallicities than galaxies with low SSFR. We thus 
derived SSFR for the 7 galaxies in our sample, using the Ho' lu- 
minosity, the SFR calibration by Argence & Lamareille ( 200^ 
[eq. 17] and the stellar mass given by the SED fitting. To be able 
to compare with the entire VVDS sample, we also derived spe- 
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Fig. 3. Specific star formation rate vs. stellar mass for VVDS 
galaxies in the 02h and 22h fields (see § l3.1l for mass estimates). 
Values are obtained with S FR derived from fOii13727A lumi- 
nosity, with the calibration of lArgence & Lamareillel(l2009l) . The 
histograms represent the mean specific SFR distribution (bold 
black line) and Icr scatter (thick black lines). The blue triangles 
stand for the seven galaxies of our sample, whi ch SFR ha ve bee n 
derived from Ha luminosity, as in lArgence & Lamareillel(l2009l) . 



cific SFR for galaxies in the 22h + 02h "wide" fields, for which 
we had a measure of F([Oii13727A) an d an estimate of the stel- 
lar mass. lArgence & Lamareill3 ( 12009 ) provide a calibration of 
the SFR based on the [On]3727A line [eq. 18] which is consis- 
tent with the Ha calibration. In Fig. [3] we show the SSFR as a 
function of stellar mass for the full VVDS 22h + 02h "wide" 
sample, together with our z ~ 1.4 VVDS galaxies. As one can 
see our galaxies lie on top of the mean distribution in SSFR, 
even out of the Icr limits (except for one), which could explain 
why their metallicity is so low. As a matter of fact our sample is 
clearly biased toward high specific star format ion rates. Several 
explanations are proposed for this trend (see Ell ison et al.l2008al 
for more details): galactic winds could be more efficient in ex- 
tracting metals at higher SSFRs, inflows of metal-poor gas could 
increase the SFR, etc. . . The fact that these seven galaxies show 
on average higher SSFRs than the global VVDS sample is not 
surprising as they were selected on the basis of the intensity of 
the [Oii]3727Aline (see Sect. 12. It . This bias could explain why 
the MZ relation we find at z ~ 1.4 seems shifted toward low 
metallicities. Another explanation, which cou ld be linked to the 
SSFR one, has been recently proposed by Peeples et al.l (l2009l) . 
These authors claim that the presence of tidal interactions in ma- 
jors mergers could possibly account for a fading in oxygen abun- 
dance as metal-poor gas may inflow from large radii. In our sam- 
ple, VVDS220596913, VVDS220544103 and VDS0201 16027 
seem indeed to undergo ti dal interactions as revealed by their 
peculiar velocity fields (see lEpinat et al.ll2009l) . 



Fig. 4. Relation between stellar mass and the N/0 abun- 
dance ratio. The solid line is the median relation of 
the SDSS star-forming galax ies in the local universe(see 
iPerez-Montero & Continill2009l fig. 14). The dashed fines indi- 
cate the scatter of this relation. The filled black squares corre- 
spond to the four VVDS galaxies for which we have measure- 
ments of the sulfur lines to derive N/0 ratio. 

4.2. Influence of the N/0 ratio 

Figure |4] shows the relation between the stellar mass and the 
nitrogen-to-oxygen abundance ratio N/0 for the SDSS star- 
forming galaxies in the local universe ( iPerez-Montero & Continil 
2009, their Fig. 13). The four VVDS galaxies at z ~ 1.4 for 
which we were able to derive \og{N/0) are also shown in this 
figure as filled squares. Even if the dispersion around the re- 
lation defined by the SDSS sample is high, it is clear that the 
four VVDS galaxies are consistent with this relation. The influ- 
ence of the nitrogen-to-oxygen abundance ratio on the metal- 
licity is sho wn on Fig. [2l with the downwa rding arrows. As 
noticed by iPerez-Montero & Continil (|2009|) the metallicities 
are significantly shifted to lower values for high N /O abun- 
dances ratios. This is especially the case for VVDS020147106 
and VVDS220014252 which exhibit the highest N/0 values. 
Looking at Fig. |2] we must keep in mind the important error 
bars going up to nearly 0.5 dex for the N/0 ratio (see Table [3]l. 
Nevertheless, even if we take into account these errors, the 
downwarding shifts in metallicity - for VVDS020147106 and 
VVDS220014252 - are too important to deny any effective de- 
cline in oxygen abundance: indeed the drops in oxygen abun- 
dance are larger than the error bars themselves, which guarantee 
an effective decline. 

5. Conclusion and prospects 

We have determined the stellar mass and metallicity of seven 
VVDS star-forming galaxies at z ~ 1.4 observed in the NIR 
with the integral-field spectrograph SINFONI. Keeping in mind 
the uncertainties inherent to the metallicity and mass estimates, 
and considering also that our small sample does not allow to 
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draw any firm conclusion, we can however draw the following 
points. Our sample of seven z ~ 1.4 galaxies shows if any, a mild 
correlation between metallicity and stellar mass (~ 0.1 logarith- 
mic slope). The mean metallicity is 12 + log(0/H) - 8.62 with a 
scatter of 0.22 dex. All the gala xies in our sample are just under 
the relation of lErb et al] (l2006h at z ~ 2, except for one galaxy 
that lies just above. The low metallicity measured in these galax- 
ies at z ~ 1 .4 could be due to their high specific star formation 
rate and/or to the tidal interactions allowing for inflow of metal- 
poor gas in the cases of VVDS220596913, VVDS220544103 
and VDS0201 16027. Accounting for the dependence of oxygen 
abundance estimate o n the N /0 abundance ratio, as suggested 
by ,Perez-M ontero & Co ntinil ((2009), would lower significantly 
the metallicity for two galaxies showing the largest N/ O abun- 
dance ratios. When completed, the MASSIV survey will gather 
a representative sample of roughly 80 galaxies observed with 
SINFONI within 1 < z < 2. It will allow to shed new light on 
these issues enabling to probe the dynamical and physical prop- 
erties of these galaxies. Concerning the MZ relation, a statistical 
study will be undertaken over a large range of stellar masses, 
computing median metallicities in bins of mass and deriving the 
median relation between 12 + log(0/H) and M* and its intrin- 
sic scatter. Comparing these observations to the p redictions of 
galaxy evolution models (e.g. Ide Rossi et al. I I2OOI will also be 
an major issue to investigate within MASSIV. 
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